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In all theoretical treatments of electron transport through single molecules between two metal 
electrodes, a clear distinction has to be made between a coherent transport regime with a strong 
coupling throughout the junction and a Coulomb blockade regime in which the molecule is only 
weakly coupled to both leads. The former case where the tunnelling barrier is considered to be 
delocalized across the system can be well described with common mean-field techniques based on 
density functional theory (DFT), while the latter case with its two distinct barriers localized at the 
interfaces usually requires a multideterminant description. There is a third scenario with just one 
barrier localized inside the molecule which we investigate here using a variety of quantum-chemical 
methods by studying partial charge shifts in biphenyl radical ions induced by an electric field at 
different angles to modulate the coupling and thereby the barrier within the 7r-system. We find 
steps rounded off at the edges in the charge versus field curves for weak and intermediate coupling, 
whose accurate description requires a correct treatment of both exchange and dynamical correlation 
effects is essential. We establish that DFT standard functionals fail to reproduce this feature, while 
a long range corrected hybrid functional fares much better, which makes it a reasonable choice for 
a proper DFT-based transport description of such single barrier systems. 

PACS numbers: 73.63.Rt, 73.20.Hb, 73.40.Gk 
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I. INTRODUCTION 

An important aspect in the vibrant new field of single- 
molecule electronics is the classification of the tun- 
nelling process with respect to the coupling between the 
molecule and metallic leads. Coherent tunnelling (CT) 
with an electron passing through the junction in one 
step (strong coupling) is contrasted with Coulomb block- 
ade (CB) with the electron resting on the molecule in a 
two-step process (weak coupling), which can be distin- 
guished experimentally^^. The theoretical description 
of electron transport through single-molecule junctions 
commonly involves the combination of a non-equilibrium 
Green's function (NEGF) formalism^ for defining the 
semi-infinite boundary conditions of the transport prob- 
lem and density functional theory (DFT) with typically 
a semi-local exchange-correlation (XC) functional for a 
reasonably accurate description of the ground-state elec- 
tronic structure of the junction at an atomistic levelA^ 
While this approach is known to reproduce the smooth 
current- voltage (I/V) behaviour characteristic of CT, it 
does not capture the distinct steps in I/V curves which 
are the signature of CB in experiments. There are two 
different schools of thought for explaining this deficiency. 
The first one relates it to the existence of spurious self- 
interaction errors (SIE) in DFT calculations with semi- 
local functionals^^— , i.e. essentialy to the approxima- 
tion for exchange in the XC functional, while the other 
blames the single determinant nature of DFT and there- 
fore its inadequacy in describing correlatio n 13 15 i 16 . 

In a way, the distinction between CT and CB made 
only from the coupling to the electrodes considers the 
molecule at the heart of the junction as a featureless 




FIG. 1: In a biphenyl molecule, the coupling between the n 
systems on its two phenyl halves is defined by the torsion 
angle between the rings: the more coplanar the stronger the 
coupling. We charge this molecule by subtracting a single 
electron and apply an external electric field for shifting partial 
charges between the two components of the resulting cation. 



quantum dot. For molecular electronics, however, it is 
particularly relevant to pay attention to the contribu- 
tion of the molecule to the transmission of the junction 
and its theoretical description. It is not unreasonable 
to consider that similar problems to those encountered 
in DFT treatments of the double barrier situation might 
also occur for electron transport through a single bar- 
rier. Therefore, we focus in this work on the description 
of a single intramolecular barrier, by studying a molec- 
ular model system and performing Hartree-Fock (HF), 
post-HF and DFT based quantum-chemical calculations 
to describe the charge shifts induced by an external elec- 
tric field. In particular, we investigate the charge shifts 
inside biphenyl radical ions, with the torsion angle be- 
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twccn the two phenyl rings acting as a parameter for the 
coupling of the two 7r-conjugated rings (see Fig. [1]). 

For the interpretation of our results we establish a link 
between our calculations and an alternative perspective 
on the limitations of DFT presented by Cohen and co- 
workers, recentl y 17 20 . In this framework DFT is con- 
sidered to be exact and universally applicable in princi- 
ple but the XC functionals currently in use are argued 
to suffer from two fundamental shortcomings, namely a 
derealization error and a static correlation erro r 17 20 . 
The former is closely linked to SIE and completely ab- 
sent in HF, while the latter affects HF even more than 
conventional DFT. For both problems, chemically sim- 
ple but methodologically still challenging test systems 
have been identified, namely the dissociation curves of 
and H2, respectively. In this picture, it follows nat- 
urally that for H^" only exchange matters because the 
system contains only a single electron, while for H2 at 
inter-atomic distances far beyond the optimal bonding 
length, the non-dynamic (or strong) correlation of the 
two unpaired electrons on the separated H atoms has to 
be correctly described. 

While Cohen et a h 17 20 and other a 21 ' 22 developed their 
ideas for the derivation of better XC functionals for gen- 
eral purposes, they also stated that their systematics 
should be relevant in particular for the description of 
electron transport problems. This suggestion has not 
been followed up by transport specific calculations so 
far while such studies appear to be necessary because: 

i) as outlined above, in electron transport theory, there 
is still no agreement within the field about the relative 
importance of exchange or correlation, respectively, and 

ii) there is still no NEGF-DFT approach available which 
can reproduce steps in I/V curves correctly in general. 
While methods exist, which go beyond a mean-field de- 
scription^—, these techniques are computationally de- 
manding and do not lend themselves easily to chemical 
interpretations in terms of molecular orbitals (MOs). It 
is therefore highly desirable to find out if there is an 
available XC functional which is able to describe steps 
in I/V curves correctly within NEGF-DFT. In a recent 
study of the Anderson junction this question is addressed 
using a Bethe ansata^l but the findings are limited to a 
single level occupied with a single electron and need to 
be complemented by highly accurate quantum chemical 
calculations on real molecular systems. We consider our 
current work to be a move in this direction, by showing 
that common density functionals fail to reproduce the 
expected intra-molecular charge shifts correctly at weak 
coupling in particular, while a long-range corrected func- 
tional behaves reasonably. 

The paper is organized as follows: in the next two sec- 
tions, we establish a direct link to the work of Cohen et 
all 17 20 by arguing that the artefacts in the description 
of charge shifts are related to the H2 and dissociation 
problems for the biphenyl neutral molecule and cation, 
respectively. In Section ITVl we present DFT calculations 
on charge shifts in biphenyl cations (we obtained quite 



similar results for anions which we do not discuss explic- 
itly in this article) using a variety of common XC func- 
tionals with the near-coupled cluster quadratic singles 
and doubles method (QCISD) as reference. We find that 
long-range corrected (LC) hybrid functional^— pro- 
vide satifactory solutions for the whole coupling range. 
Finally, we provide a summary where we also discuss the 
relevance of our results for NEGF-DFT-based electron 
transport calculations. 



II. SHIFTING CHARGES INSIDE MOLECULAR 
CATIONS AND NEUTRAL MOLECULES AS A 
PROBE FOR COUPLING IN MODEL SYSTEMS 

In order to establish a link between the dissociation of 
molecules and ions on the one hand and electron trans- 
port through single molecule junctions on the other hand, 
we take an intermediate step and pose an associated 
question, namely how dissociation can be linked to shift- 
ing charges inside single molecules with an electric field. 
First we discuss how such charge shifts should evolve in 
molecules or ions composed of two identical halves, such 
as in H2 and biphenyl (BP). If no field is applied, the 
MOs contain equal contributions from atomic orbitals 
(AOs or, in general, fragment orbitals which we will also 
call AOs for the sake of simplicity in the following) of 
both parts by symmetry, and hence the charge is equally 
distributed regardless of the strength of the coupling. If 
a longitudinal electric field (EF) is introduced the left 
and the right sides experience different potentials, the 
AO levels are only shifted upwards in energy on one side 
and downwards on the other side, which will influence 
the composition of the MOs accordingly (in a first ap- 
proximation neglecting the polarization of AOs). 

When the coupling between the parts is sufficiently 
strong, the impact of the applied field will be an unequal 
contribution of the left and right AOs to all MOs, that 
is a polarization of the system, so that the composition 
of the MOs will change smoothly with EF. This grad- 
ual polarization of occupied MOs, in particular of the 
HOMO, induces a continuous shift of fractions of an elec- 
tron from one side to the other. The weak coupling limit 
is more challenging for the description of charge shifts 
inside molecules. For zero coupling between the parts, 
in the absence of EF, all the MOs of our two-component 
system come in degenerate pairs arising from a symmet- 
ric and antisymmetric linear combination of the left and 
right AOs of the same type and also equal in energy to the 
respective AOs on both parts, since hybridization effects 
are ruled out by definition. If the symmetry constraints 
arc lifted by applying EF, the MOs in every pair are split 
and take the shape of their constituent AOs. In the weak 
coupling limit there are fundamental differences between 
systems with even and odd numbers of electrons, which 
we discuss below with the help of Figs. [5] and [3] 

For the neutral H2 and BP molecules, which contain an 
even number of electrons, the process of shifting charges 
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involves mainly the HOMO and the LUMO which in the 
zero coupling limit are reduced to the corresponding AOs 
on each half. Depending on the intensity of the applied 
field, there are two possible ground states obtained either 
by populating each orbital with one electron, or by popu- 
lating the lower one with both electrons. For H2 (Fig. [2]), 
the open-shell configuration with a field applied (panel E) 
has the same level occupation as the zero-field ground- 
state (panel B) with one electron on each side until the 
potential difference AV (induced by EF and increasing 
with it) between the AO levels exceeds the charging en- 
ergy AU. At this point the system is transformed to a 
closed shell configuration with both electrons on one side 
(panel D), where an electron "jumps" from the higher ly- 
ing singly occupied level to the lower lying one, making 
the latter doubly occupied. 

We use a model system with four electrons in the fron- 
tier orbitals to mimic the behaviour of the 7r system of 
a neutral BP molecule with the phenyl rings perpendic- 
ular to each other as indicated in panels B, D and E of 
Fig. H where only the local HOMOs and LUMOs of the 
phenyls are taken into account). The zero- field ground 
state (panel B) is closed-shell, and EF (shifting all the 
occupied and vacant levels at each side in the same di- 
rection) has to reduce the gap (A gap ) between the right 
HOMO and the left LUMO (panel E) until an electron 
transfer between these two levels creating two SOMOs 
at the opposite sides (panel D) becomes favorable. In 
this case, the electron jump leads to a transition from 
a closed to an open shell configuration. We note that 
the relation between AU and A gap with respect to AV 
depends on the theoretical framework in this case, since 
the charging energy is contained in the energy gap for 
HF but contrarily it needs to be explicitly added to the 
Kohn Sham gap for arriving at addition energies within 
DF T 32 ' 33 . For our purpose in this article, this particu- 
lar distinction is not relevant since we focus on cations. 
In order to summarize the common features between the 
two cases with an even number of electrons, we conclude 
that for the electron to jump from one side to the other 
a charging energy AU has to be overcome by the applied 
field EF. 

For the cations H^ and BP + with one and three elec- 
trons in the model system, respectively, for zero coupling 
and zero field, a single electron occupies one MO of the 
degenerate pair (see the upper left side panels of Figs. [H 
and [3]) . Due to symmetry, the probability to find the 
electron on either of the two AOs (halves of the cation) 
is equal. By applying EF this symmetry is lifted along 
with the degeneracy of the AOs and the electron becomes 
localized on one side. The spatial charge distribution 
then crucially depends on how localized or delocalized 
the fragment orbitals are. It is essential for the argu- 
ments in our paper that for the breaking of symmetry 
caused by the physical boundary conditions imposed by 
the external field an infinitesimally small AV is suffi- 
cient and no charging energy AU has to be provided or 
in other words in the zero coupling limit no resonance 
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FIG. 2: Schematic orbital diagram of the shift of electrons in 
Hj" (panels A and C) and H2 (panels B, D, and E) for the limit 
of zero coupling between the orbitals before (A, B) and after 
(C, D and E) the application of a longitudinal electric field. 
For H2 a critical voltage AV exceeding the charging energy 
AU has to be applied (D), while an infinitesimally small AV 
is enough for H^~ (C). The corresponding MO patterns are 
shown below their energy levels, where positive and negative 
signs of the wavefunction are illustrated as blue and white 
spheres, respectively, and the position of the nuclei marked 
by black dots. 
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FIG. 3: Schematic orbital diagram of the shift of electrons 
in a longitudinal electric field between the ir systems of the 
two phenyl rings in a neutral biphenyl molecule (panels B, 
D, and E) and cation (panels A and C) for the limit of zero 
coupling. The portrayed orbitals are the local HOMOs and 
LUMOs on both rings, respectively. While for the neutral 
molecule the threshold AV has to be exceeded (incorporating 
both the HOMO-LUMO gap and the charging energy) for a 
full charge transfer (D), an infinitesimally small electric field 
is enough for the cation (C). 
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energy is lost by localization in the more stable of two 
quasi-degenerate structures. This distinguishes the sys- 
tems with an odd number of electrons from those with 
an even number significantly. 

In order to elucidate the connection between the two 
different physical phenomena of shifting charges inside 
molecules or ions and their dissociation by stretching 
bonds, it is worth noting that there are two key param- 
eters in both cases determining the impact of their elec- 
tronic structure: (i) the relative position of the AO levels 
of the isolated fragments and (ii) the coupling between 
their AOs. In the course of dissociation, for a chosen 
molecule, (i) remains constant being determined by the 
chemical nature of the fragments, while (ii) is continu- 
ously varied being a function of the distance. The pa- 
rameter landscape of (i) can also be scanned in the case 
of dissociation, although discontinuously, by comparing 
the curves of various molecules. In contrast, when EF 
is applied to a given molecular system, (ii) remains con- 
stant (in a first approximation, neglecting the polarizing 
effect of the field on the shape of AOs), while (i) is var- 
ied continuously. If this is done for various couplings, it 
means that the scope of (ii) is also probed. Therefore, the 
entire parameter space explored in the description of the 
electronic structure of dissociating molecules is equiva- 
lent to that when shifting charges inside molecules with 
a longitudinal electric field and a variable torsion angle 
to modulate the electronic coupling. It is therefore rea- 
sonable to assume that methodological findings on the 
dissociation of ITj~ and H2 are also relevant for the ap- 
plicability of the investigated methods to the problem of 
charge shifting with an external field. This conjecture, 
however, has to the best knowledge of ours never been 
explicitly tested which is a key motivation for the work 
presented in this article. 



III. FROM MODEL SYSTEMS TOWARDS DFT 
CALCULATIONS 

When performing DFT calculations of the dissociation 
potential energy curves of H2 and H^~ , respectively, two 
distinctly different scenarios can be identified with regard 
to the total energy for a large separation of the two H 
atomsi£. 

For Hj, DFT finds the state with the single electron 
in the system delocalized on the two H positions al- 
ways overstabilizcd (derealization error) vs. the state 
in which it is localized at one of the hydrogens, while 
in the low coupling (dissociation) limit these situations 
should tend to degeneracy. This deficiency results in dis- 
sociation to a state with fractional charges (which in itself 
is not an error for H^, but can lead to errors for other 
diatomics^i) and a binding energy that is far too low^i. 
It has been shown that for obtaining the correct physical 
relation between the total energies of atoms with integer 
and fractional charges, the XC-potential would need to 
exhibit a derivative discontinuity, whereas the common 



semi-local functional do not^— . Since this problem 
is related to SIE, HF calculations do not suffer from it 
(the self-interaction contributions of the Hartree and ex- 
change potentials cancel out exactly), and therefore the 
binding energy of H^ can be obtained exactly with HF. 

For H2 the situation is very different. Here, the prob- 
lem for mean-field computations is due to an insufficient 
description of the correlation between the spins of the 
two electrons^ in both DFT and HF (static correlation 
error). The latter case which corresponds to shifting 
charges in a neutral molecule (see the right side pan- 
els of Figs. [2] and [3]) closely resembles the electron jumps 
between molecules which we described in Ref 

As our aim is to model the electron transport through a 
single intra-molecular barrier, we argue that the cations 
on the left sides of Figs. [2] and [3] are closer to such a 
setup than the neutral molecules, because an electron 
moving through a junction is provided by the electrodes. 
The injection of a carrier into the molecule can involve a 
charging energy but only if the molecule is weakly cou- 
pled to the electrodes. This is not the case we want to 
investigate in our present article, where we are only con- 
cerned with the effect of a single intra-molecular barrier 
and therefore assume that the electron (or in our study 
the hole, which is equivalent) has already reached the 
molecule and hence focus on cations. 

For H^ and BP + , which are open-shell systems as 
a whole, with and without the electric field, a closed- 
shell description becomes formally inapplicable. How- 
ever, while the correlation effects are strictly inexistcnt 
for H^ , which allows unrestricted HF (UHF) to describe 
it properl y 21 i 38 , the situation can be more complicated 
for BP + . Although static correlation which is the main 
source of error for H2 (and neutral BP) should be absent 
for the biphenyl cation as well as for H^ , BP + is not a 
one-electron system, so that dynamical correlation can 
become significant due to the interaction of the electron 
in the singly occupied orbital with all other electrons in 
the radical ion. Such effects have already been identi- 
fied in the simpler but related case of allyl radicals^ and 
have also been discussed from the perspective of artifi- 
cial symmetry breaking or the competition between spin 
polarisation and spin derealization^. 

The question now is how DFT with different XC- 
functionals deals with this situation. It is not easy to 
predict a priori and addressing this issue is the main ob- 
jective of this article. Because the derealization error is 
known to artificially stabilize the delocalized system for 
dissociation, leading e.g. to spurious fractional charges 
on the fragments^, it is natural to assume that for EF in- 
duced charge shifts, the same derealization error makes 
the polarisation of a molecular cation more gradual (or 
in other words the charge distribution between its halves 
more equal) , since in both cases the deficiency is due to an 
overcstimation of the coupling throughout the molecule. 

For the first-principles calculations presented in the fol- 
lowing, H^ is an impractical system because the coupling 
between the atoms, which is the crucial quantity to probe 



5 



the connectivity of a frontier MO relevant for electron 
transport, can only be varied by their distance. For BP + 
on the other hand, the coupling between the local ir elec- 
trons can be controlled by varying the torsion angle be- 
tween the rings, which can also be adjusted by chemical 
means thereby allowing for its use as an experimentally 
trackable coupling parameter for electron transport mea- 
surements in single molecule junction o 41 ! 42 . 



IV. FINDING A SUITABLE XC-FUNCTIONAL 
FOR SHIFTING CHARGES IN A BIPHENYL 
CATION 

All numerical results presented in this section are 
based on first-principles calculations performed with the 
GAUSSIAN09 package^ 3 -. We use this code because it 
allows for a direct comparison (i.e. within the same 
computational setup and using the same basis sets) of 
DFT based techniques with semi-local or hybrid XC- 
functionals with HF and post-HF implementations of 
wavefunction theory. In Fig. we plot charge vs. 
voltage curves calculated with UHF (black lines) for a 
biphenyl cation for a variety of torsion angles ranging 
from 110° to 170° in the optimized geometry of a neu- 
tral BP molecule with the respective torsion imposed as 
a constraint, i.e., from very weak to very strong cou- 
pling between the n systems on the two phenyl rings, 
respectively, as obtained from a summation of Mullikcn 
charges. We note that for weak coupling (110°-120°) not 
only the smallest field we apply already breaks the sym- 
metry between the two components and the charge lo- 
calizes on one of the two phenyl rings but the zero-field 
charge distribution also has a broken symmetry, because 
the symmetric UHF solution turns out to be unstable, 
i.e. higher in total energy than the asymmetric one. The 
singly filled MO (SOMO) in particular always tends to 
localize preferentially on one of the two rings. For strong 
coupling (170°), this SOMO is well delocalized over both 
components also for high fields and charge is therefore 
transferred continuously by a polarisation of this orbital. 
As can be seen in Fig. Hk, this polarizability gradually 
increases as the driving force derealization decreases for 
higher coupling/torsion angle. 

Although UHF describes exchange interactions ex- 
actly, its symmetry breaking problem at zero field makes 
it a questionable bechmark even for calculating partial 
charges at a finite voltage where the field is breaking 
the symmetry as a physical boundary condition anyway. 
Though counterintuitive, this is due to the continuity of 
the finite field solution with the artificially asymmetric 
zero field solution. This implies that if the zero field so- 
lution is wrong, the finite field solution has to be wrong 
as well. At zero field and with the positions of the nuclei 
fixed the geometry of the wavefunction has to be sym- 
metric because of the symmetry of the molecule, so that 
we need to find a method providing such a solution as a 
benchmark. Sometimes this can be achieved by increas- 



ing the size of the basis set but we checked that it is 
not the case here. Since the complete active space self- 
consistent field (CASSCF) and configuration interaction 
singles and doubles (CISD) methods still break electronic 
symmetry without bias to a significant degree, we chose 
QCISD calculations as a benchmark, a method for which 
it has been demonstrated for allyl radicals that higher or- 
der correlation effects are well described 4 ^. We checked 
that for our system the symmetry breaking tendency of 
this method is negligible and found that its main effect 
on the charging curves compared to UHF at an angle of 
110° is to round off the edges of the steps considerably. 

DFT calculations within the local spin density approx- 
imation (LSD) and using the Vosko-Wilk-Nusair (VWN) 
parameterisation 4 ^ for the XC-funtional do not suffer 
from the symmetry breaking problem but its deviation 
from the physically correct curves are equally serious. 
Within this framework the partial charge curves become 
linear regardless of the torsion angle (red lines in Fig.^J,), 
which is due to the lack of a derivative discontinuity in 
the XC-potential, and leads to an overstabilisation of 
fractional charges on both phenyl rings compared to a 
distribution of integer charge o 18 i 35 . This deficiency can 
also be understood in terms of a SIE-induced derealiza- 
tion erro r 45 ' 46 , even for a torsion angle of 110°, where 
the SOMO is almost equally ditributed on both rings at 
high fields with LSD- VWN, since delocalized states are 
always lowered in energy artificially. 

In Fig. UJd wc compare charging curves for an indi- 
vidual phenyl ring within a biphenyl cation at a tor- 
sion angle of 110°, as calculated with a variety of dif- 
ferent XC-functionals and measure them against QCISD 
as a physically correct benchmark. It can be seen that 
the generalized gradient approximation (GGA) with a 
Perdew-Burke-Ernzerhof (PBE)i£ parameterisation for 
the XC-functional as well as the hybrid exchange ap- 
proach B3LY P 48 ' 49 provide only incremental quantitative 
improvements when compared with LSD- VWN. BHand- 
LYP, with 50% of HF exchange, approaches the QCISD 
curve somewhat closer, still underestimating significantly 
the polarization induced by the field. We find, how- 
ever, a quite different behaviour with the long-range 
corrected hybrid functional LC-a;PBE 29 31 which, how- 
ever, somewhat overshoots with respect to the QCISD 
polarisation at 110°. As a further more stringent test, 
we check whether a LC-wPBE parameterisation of the 
XC-functional reproduces the QCISD results also for a 
more intermediate regime. For this purpose, we show 
charging curves for torsion angles of 110° and 130° in 
Fig. 0b, where qualitatively we find that DFT-LC-wPBE 
correctly reproduces the gradual change in curve shape 
from rounded steps to continuously linear with an in- 
crease in coupling between the ir electrons of the phenyl 
rings, while BHandHLYP does not capture this progres- 
sion and exhibits straight lines for both angles. 

In LC-wPBE, the exchange component of the Hamil- 
tonian is separated into a short-range (SR) and a long- 
range (LR) contribution with a range-separation param- 
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FIG. 4: Partial positive charges on one of the two phenyl rings in a biphenyl cation induced by an electric field as calculated: 
a) for various torsion angles between the rings with UHF (black lines) and DFT making the LSD approximation in a VWN 
parameterisation (red lines); b) at a torsion angle of 110° with UHF (black), QCISD (magenta) and DFT with a variety of 
different XC-functionals, namely LSD- VWN (red), GGA-PBE (orange), B3LYP (green), BHandHLYP (cyan) and LC-^PBE 
(blue); c) for torsion angles 110° and 130° with QCISD (magenta) and DFT using LC-wPBE (blue) and BHandHLYP (cyan) 
as parameterisation of the XC-functional. 



eter uj which was empirically fitted to the Gaussian de- 
fault value of 0.4 Bohr -1 (Ref.— ), which can be further 
refine d 50 ' 51 . This scheme, which uses pure PBE for the 
SR part and 100% HF exchange for the LR part, was 
originally developed to correct the long-range behaviour 
of the XC-potential for molecules^ but has also been 
recently found to reduce SIE systematically in systems 
with fractional electron numbers 3 -!; following the logic of 
Cohen et alJ&, this should also improve the description of 
the dissociation curve of molecular ions such as H^. The 
explanation for this benign property of the functional 
offered by the developers of LC-wPBE is that while a hy- 
brid functional with less than 100% of HF exchange does 
not satisfy any universal constraint beyond those satisfied 
by the underlying semi-local functional, the LR-corrected 
LC-wPBE in contrast recovers the exact asymptote of the 
exchange potential in molecules, which is expected to im- 
prove the description of the density tail regions 3 -!. It is 
exactly this LR asymptotic behavior of the exchange po- 
tential which matters for both H^~ dissociation and the 
relative stability of fractional charge d 35 ' 36 . Our study 
thus indicates that long-range corrected hybrid function- 
al such as LC-wPBE are best suited for investigating the 
electron transfer through biphenyl molecules with near 
perpendicular orientation of the rings. To the best of our 
knowledge, such functionals although proposed for other 
purposes have not yet been implemented in any of the 
common NEGF-DFT codes. 



V. SUMMARY AND DISCUSSION 

In summary, we established a connection between 
intra-molecular electron transfer in biphenyl cations and 
neutral molecules and the well-studied dissociation prob- 
lems for H^" and H2. It is known for the former that 
DFT methods with semi-local XC functionals suffer from 
a derealization error - a deficiency of the exchange part 
of the Hamiltonian - while for the latter case, they are 
inadequately describing static correlation effects. In or- 
der to make such a link, we took an intermediate step 
and considered first schematically the electron transfer 
inside H2 and H^~ as well as in the biphenyl molecule 
and cation with a weak coupling between the compo- 
nents. Thereby, we could make a distinction between the 
case: (i) of H2 and the neutral biphenyl molecule, where 
charges have to be created locally in an originally closed 
shell system and static correlation errors dominate (we 
have already studied this case extensively in our previ- 
ous worki£) and (ii) H^ and the biphenyl cation where 
an already existing charge is just localized on one com- 
ponent by symmetry breaking upon application of the 
electric field and exchange; for BP + also dynamical cor- 
relation needs to be correctly described in order to avoid 
the derealization error and artcfactual discontinuities in 
the charging curves. In order to further investigate the 
latter scenario, we performed quantum-chemical calcu- 
lations on a biphenyl cation with various torsion angles 
and showed that the long-range corrected hybrid func- 
tional LC-wPBE is a good parameterisation for the XC- 
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functional within DFT in this context. 

In this contribution, we investigated the transfer of 
partial charges from one component of a molecular cation 
to another, as induced by an external electric field. Al- 
though such a study gives an insight into the ability of 
a theoretical approach to describe the transparency of 
a coupling barrier inside a molecule which is relevant 
for electron transport through a single molecule junc- 
tion, it is clearly not directly comparable to a theoretical 
description of the transport problem as provided by a 
NEGF-DFT scheme^, where the electrodes are treated 
explicitly as metal slabs made semi-infinite by Bloch-type 
conditions. For the situation of a single intra-molecular 
barrier discussed in our article, however, the similarities 
between intra-molecular electron transfer and electron 
transport through a junction in terms of methodologi- 
cal requirements are significant. In the case of charge 
transfer inside a biphenyl cation, no local charges are 
created by the field and therefore no charging energy is 
involved. Similarly, in a transport setup for a biphenyl 
molecule bonded to gold electrodes by thiol anchors^- 
— , assuming that the thiol groups ensure ideal coupling 
of the biphenyl MOs to the leads, perpendicular rings 
would create a situation where two semi-infinite electron 
densities are separated by one coupling barrier. In such a 
system there would be no electron localisation on one of 
its components with a finite extension and therefore no 
charging energy U, which suggests that an approach with 
an accurate description of exchange and dynamical corre- 
lation such as NEGF-DFT with LC-wPBE for which we 
make the case in our article but also the GW method^ 
might be sufficient for a proper treatment of this sce- 
nario. Interestingly, we also predict that a method which 
is good at capturing static correlation such as CASSCF 
fails to describe dynamical correlation correctly in such 
a single barrier setup. 

As validation for our assumption we refer to Ref^, 
where the transmission of biphenyl with thiolate and car- 
bodithiolatc anchors (BDT and BDCT, respectively) was 
calculated as a function of the torsion angle with NEGF- 
DFT. In this work, it was found that the zero-bias con- 



ductance with two coplanar phenyl rings differs almost 
by an order of magnitude with BDCT being more con- 
ductive. However, as the torsion angle approaches 90°, 
not only the conductance of both molecules in the junc- 
tion decreases, but also the difference between BDCT 
and BDT fades out in such a way that the conductance 
becomes practically the same for the perpendicular con- 
formation. We would interpret this finding as an indi- 
cation that, as the coupling between the 7r-systems of 
the phenylenes decreases, the limiting barrier for con- 
ductance gradually shifts from that determined by the 
anchors (and favoring BDCT) to that determined in- 
tramolecularly by phenylenes (and thus equal for both). 
This vindicates our assumption above that in a trans- 
port setup for such a scenario the anchors to the leads 
can be regarded as ideal contacts, which makes it in effect 
a single-barrier system. This observation can serve as a 
justification for concluding from our molecule-centered 
approach to charge transfer that in electron transport 
calculations with NEGF-DFT rounded off steps would be 
found in the I/V curve for a single molecule junction with 
a substantial intra-molecular barrier (weak coupling) if 
an LC-wPBE (or equivalent) functional is used to ensure 
a proper exchange description. In order to observe this 
feature experimentally, a three-electrode scenario with a 
gate in addition to source and drain would be required 
(which is the standard for measurements on single elec- 
tron transistors) for justifying our assumption that the 
molecule is already charged when the tunnelling process 
through the single barrier is initiated by the field. 
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